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We discuss the change of the Kondo effect in the Josephson junction formed by the indirect cou¬ 
pling between a one-dimensional Dill-class topological and s-wave superconductors via a quantum 
dot. By performing the Schrieffer-Wolff transformation, we find that the single-electron occupation 
in the quantum dot induces various correlation modes, such as the Kondo and singlet-triplet corre¬ 
lations between the quantum dot and the s-wave superconductor and the spin exchange correlation 
between the dot and Majorana doublet. Moreover, it plays a nontrivial role in modifying the Joseph¬ 
son effect, leading to the occurrence of anisotropic and high-order Kondo correlation. In addition, 
due to the quantum dot in the Kondo regime, extra spin exchange correlations contribute to the 
Josephson effect as well. Nevertheless, if the Dlll-class topological superconductor degenerates into 
D-class because of the destruction of time-reversal invariance, all such terms will disappear com¬ 
pletely. We believe that this work shows the fundamental difference between the D- and Dlll-class 
topological superconductors. 

PACS numbers: 73.21.-b, 74.78.Na, 73.63.-b, 03.67.Lx 


I. INTRODUCTION 

Observation of the Kondo effect in semiconducting quantum dots (QDs) has opened a new area of research 
about the quantum transport through nanoscale and mesoscale systems»i Since QDs have advantages to 
couple and then form QD molecules, the Kondo effect in multi-QD systems exhibits various forms. For 
instance, the spin and orbital Kondo effects, two-stage Kondo effect, and Kondo-Fano effect have been 
observed, which induce abundant transport behaviors.— In addition, the Kondo effect plays a special role in 
affecting the Andreev reflection in the systems with superconducting leads. In one Josephson junction with 
an embedded QD, the interplay between the Kondo and Josephson effects can efficiently drive the Josephson 
phase transitionj^i^ Namely, if the Kondo temperature Tk is much greater than the superconducting gap 
As, the Kondo screening will dominate the system and the Josephson current will be at its 0 phase. Instead, 
if Tk ^ As, the ground state is a BGS-like singlet state, so the 7r-junction behavior occurs. 

In recent years, topological superconductor (TS) has become one main concern in the field of mesoscopic 
physics because Majorana modes appear at the ends of the one-dimensional TS which are of potential appli¬ 
cation in topological quantum computation/-^!’ Owning to the presence of Majorana zero mode, fractional 
Josephson effect comes into being, manifested as the 47r-periodic oscillation and parity-related direction of 
the supercurrent. As a typical case, in the D///-class TS, which is time reversal invariant, Majorana mode 
appears in pairs due to Kramers’s theorem?ii“— Hence at each end of the Dlll-class TS nanowire there 
will exist one Majorana doublet— Since the Majorana doublet is protected by the time-reversal symme¬ 
try, the period of its driving Josephson current becomes related to the fermion parity of the Dlll-class TS 
junction,— What’s more interesting is that in the junction formed by the coupling between the D///-class 
TS and s-wave superconductor, Josephson effect presents apparent time-reversal anomaly phenomenon.^- 

In view of the special property of the Dlll-class TS, we anticipate that in the its-existing Josephson 
junction, the Kondo effect can arouse new physics results. With such an idea, in this work we aim to 
investigate the new feature of the Kondo effect in a hybrid Josephson junction, i.e., the junction formed by 
the indirect coupling between a Dill-class TS and s-wave superconductor via one QD. Our investigation 
shows that the single-electron occupation in the QD induces multiple correlation modes, including the 
Kondo and singlet-triplet correlations between the QD and the s-wave superconductor and spin correlation 
between the QD and Majorana doublet. Besides, the Kondo QD modifies the Andreev reflection between 
the superconductors, leading to the occurrence of anisotropic and high-order Kondo correlation and extra 
spin exchange correlations which contribute to the Josephson effect. 

The rest of the paper is organized as follows. In Sec. II, the model Hamiltonian to describe electron 
behavior in the hybrid junction is introduced first. Then its low-energy effective form is derived. Then, 
different correlation modes in this structure are presented. In Sec. HI, the Schrieffer-Wolff transformation is 
carried out. In Sec. IV, the Josephson-Kondo effect is discussed. Finally, the main results are summarized 
in Sec. V. 
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FIG. 1: (a) Schematic of a QD-embedded Josephson junction with Z)///-class topological and s-wave superconduc¬ 

tors. (b) Illustration of the quantum states this junction. 


II. THEORETICAL MODEL 

The considered Josephson junction is illustrated in FiglHa), where a one-dimensional Dlll-dass TS couples 
to one s-wave superconductor via one QD. The system Hamiltonian can be written as H = Hp+Hg+Hu+HT- 
The first two terms denote the Hamiltonians of the DIII-c\ass TS and s-wave superconductor, which can be 
expressed as^ 

Hp = —/Tp cj^Cjcr — t + h-C.) -f —ICTl (72 -|-h.C.], 

j(J j<T jaa' 

Hs = ( 1 ) 

k(7 k 

Cj^ and fl^ {cja and fka) are the electron creation (annihilation) operators in the Dill-class TS and s-wave 
superconductor, respectively, with the spin index tr. ai {I = 1,2,3) is the Pauli matrix, /ip and t denote 
the chemical potential and intersite coupling in the Dlll-class TS, and ^k is the electron energy at state 
\ka) in the s-wave superconductor. Ap and As are the Copper-pair hopping terms in the two kinds of 
superconductors. Hu describes the Hamiltonian of the QD. As a typical case, we consider one single-level 
QD to be embedded in the junction and then 


Hu = eo'^Ua + U n-t-nj, (2) 

(7 

with Tier = d^da- dj. and da are the creation and annihilation operators in the QD, Eq is the QD level, and U 
denotes the intradot Coulomb repulsion. Ht describes the couplings between the QD and superconductors. 
Its expression can be given by 


Ht = -St E + E + hx., (3) 

(7 k<7 

where St and 14 are the QD-superconductor coupling coefficients, respectively. 

For the case of one infinitely long DHI-dass TS, one Majorana doublet will form at its end, we can 
therefore project Hp onto the zero-energy subspace of Hp. As a result, H can be simplified as 

H = Y,e,dlda + Un^n^ + Y.^kflfka + {Y.l^sfyUi + Y.^kflda 

<7 k(7 k k(7 

-'^Ste'-ijoada+d.c.). (4) 

<7 

joa is the Majorana operator, which obeys the anti-commutation relationship of { 700 -) 7oo-'} = ‘^Saa’- Based 
on such a form of H, the distribution of relevant quantum states in this system can be illustrated in Fig. 1(b). 

It is known that Kondo effect arises from the antiferromagnetic correlation between the localized and con¬ 
duction electrons. In an Anderson-typed system, the Kondo physics can be well described by transforming it 
into one s-d exchange model by performing the Schrieffer-Wolff transformation. For our considered junction, 
we would like to discuss the Kondo effect using such a method. The detailed discussion is shown as follows. 
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III. SCHRIEFFER-WOLFF CANONICAL TRANSFORMATION 


Since the presence of a s-wave superconductor, it is necessary to introduce the Bogoliubov unitary trans¬ 


formation to diagonalize its Hamiltonian by defining 
Hamiltonian can be written into H = Hq + Ht, where 


/fcf 


Uk 

* 

1 



. -f-ki . 


_ Vk 

* 




. And then, the system 


Ho = Sddldcr + Un^ni + Eka\^aka, 

<y kcr 

Ht = Vfc('Ufc(i],afco-+ h.c.) - -I--I-h.c.) - + h.c.) (5) 


k(7 


with Ek = It is readily found that Ho is the diagonalized Hamiltonian including the QD and 

s-wave superconductor, whereas Ht describes the hopping between the QD and superconductors. In order to 
perform the Schrieffer-Wolff canonical transformation, we need to project the Hilbert space of this junction 
into two sectors, namely, a low energy subspace and a high energy subspace, with projection operators 
Pl = \L){L\ and Ph = \H){H\ = 1 — Pl- Following the definition of projection operators, the total 
Hamiltonian can therefore be divided into the diagonal and off-diagonal parts 


Ho 


Ht 


PlHPl + PhHPh 


PlHPh + PhHPl 


Hl 0 
0 Hh\^ 

0 vf 

V 0 ■ 


( 6 ) 


Next, via a canonical transformation, H can be rotated into a block-diagonal form, i.e., 

H. According to the Baker-Campbell-Hausdorff formula, the transformed Hamiltonian H 
into series 


Hl 

AV 

can 


.-s _ 


AVf 
Hh 
be expanded 


H = H+[S,H] + ^[S, [S, H]] + [5, [^, i/]]] + • • • . 


( 7 ) 


We choose the generator S of the canonical transformation such that [5, Ho] 
part Ht is eliminated in lowest-order term. And then. 


H 


Ho + ^[' 5 ', Ht] + y]] 


n 


n=2 


(n -b 1)! 


[S,HT]r 


—Ht, hence the off-diagonal 


( 8 ) 


This result indicates that the tunneling processes can be described by second(and higher)-order terms. 
Within the low-order approximation, the effective Hamiltonian shows the form as H = Ho + AH with 
AH = A [S', Ht]- It should be understood that AH contributes to the occurrence of the Hondo effect. 

In what follows, we aim to solve the analytical form of AH by deriving the generator S. S should be one 
anti-hermitian generator obeying the relationship that — S = S^, it can thus be written into the matrix form 

. Substituting this form into equation [S, Ho] + Ht = 0, we can obtain the result 


0 -St 
S 0 


S = 


E 

\L)-\H) 


Ht 

Hh — Hl 


Pl 


h.c.]. 


(9) 


For a single-level QD, it possesses four electron states, i.e., |0), ]na), and It is suitable to introduce 

a set of projection operators in terms of occupation configuration on the QD, which are Po = (I — n-|-)(I — n^,), 
Pl = n-|-(l — nj,) -b n 4 .(l — n^), and P 2 = n-i-nj.. In the case of < 0 and Ed -\-U >0, the Hondo effect 
have an opportunity to come into being. In such a case, the ground state of the Anderson-typed system 
is a local moment jl) = ]na) configuration, whereas the high-energy intermediate states corresponds to jO) 
and [n-i-nj.). As a result, the projector into the low(high) energy subspace is Pl = Pi {Ph = Pq + P 2 )- 
In view of the expression of Ht, one can find that three kinds of coupling manners occur between the QD 
and superconductors, i.e., Vkiukd^aka + h.c.), H^^ = - J2k ^k{vl.d{^al^ + + h.c.), and 

= —St^^{e^2joadcr -bh.c.). is traditional mixing term between the conduction quasielectrons and 

d-electrons. It shows that when a conduction electron or hole is excited into the localized d-state to create 
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these excited state configurations, the corresponding excitation state is either a state with one conduction 
quasi-electron with its energy or a state with two d-electrons and one conduction quasi-hole. In these two 
processes, the excitation energies are Ai?io = Ej. — Sd and Ai?i 2 = Sd + U — E^, respectively. Accordingly, 
can be given by 


H,LM 

= VkUk{ 

k(7 


Hh-Hl 


l-ng 


Afc £d E}^ Ed U 


)al„dg-h.c.. 


( 10 ) 


Hdjd describes the scattering process between d-electrons and conduction quasiholes and vice versa. In this 
processes, the virtual excitation energies are AAqi = —Ek — Ed and AA 12 = £d + U + Ek, and then can 
be written as 

oil ^ T/ro K4°^ii + '<Jkakid^){n^-ni)'^ {vkUk^di + vldlal^){n^ - nif 

O = — 2_^ Vki-rH -- 77 - 77 -— il-C.J 


H,L,k<T 




Hh-Hl 

Ti^ 


Hh-Hl 


k<7 


Ek + Ed Ek £d-\-U 


^E^kadg h.C.. 


( 11 ) 


Next, results from the mixing between Majorana zero mode and the QD. It contributes to the virtual 
excitation energies AAgi = —Ed and AE 12 = £d + U. Consequently, can be expressed as 

gin 




H.L,< 


Hh-Hl 


= -I-_ h.C.. 


-Sd 


-£d - U 


( 12 ) 


Up to now, the generator of the Schrieffer-Wolff transformation in such a hybrid structure has been solved, 
i.e., 


kcT ka 

n. 


Ek -f Ed Ek Ed EU 


^UkadfT 


E C i — / 1 — He 

5te^{ -^ -h 


-Ed -Ed - U 


)7o<TdtT — h.C.. 


(13) 


It is readily verified that S obeys the necessary conditions = —S and [S', i?o] + Ht = 0. 

Following the expression of S, the analytical form of AH can be derived using the relation of AH = 
1[S, TSt]- Note that since the Kondo effect occurs under the condition of < 0 and Ed + U > 0, it is 
necessary to calculate AHll, which is given by AHll = ^Pl[S,Ht]Pl- After a simple derivation, one can 
find that 


AHll = ^ PlHt „ HtPl = = hI 


\H) 


Hl-Hh 


El — Eo 


El — E 2 


l^Hoi+Hl,^H2i, (14) 


APio 


AE 12 


due to Pl = Pi and Ph = Pq + P 2 - According to the result in Ea. dT^ . we can readily obtain the detailed 
form of each quantity in this equation. To be specific, H^ = 'Lka ^kUkOL\^dn{l — rig), AE(q = Ek — £d', 
Hqi = -VkVkC(kadg{l - ng), AE{1 = -Ek - Ed] HqL = Ste’-^jogdg{l - ug), AE{q^ = -Eg. Meanwhile, 
^21 = J2ka VkUkdlugUka, AEI 2 = Ed + u - Ek] H^l = - VfcU^dl , AUf^ = Ek + £d + U] H^L = 
6te~^^d^ngEoa, AE{{^ = —Ed — U. Substituting these quantities in Ea. dT^ . we can readily find that each 
term of AHll consists of four parts. To be concrete, the four parts of its first term are 


(!) : E 

kk' ^aa' 


2 

Uk ^ U 


^k ^k' 


akcrCel,, -I- 


i „t „ 

■ H-- I o^kdO^k'd' 


^k ^d ^k' ^d 


6 ^f^dlda + - + Sgn(cr)5 


(77): - y] VkVy 


UkVk UyVk' 

Ek — Ed -Ek’ — Ed 


O-kdO-k'd' + 


UkVk UyVy \ t E 

-Ek -Ed Ey -Ed / 


Syddidd + Q -I- Sgn((T)S'' 
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(III) : ^ Vfe 


f UkSte^^ ^ Uk&te^'' 

yEk-ed -Ed 


^ka'yOa' 


VkSte 


— Ek — Ed 


*•’ vl&tE'- 
-Ed 


aLlOa' 


Uk&te *2 UkSte *2 


-Ed 




Ek — Ed 

S^,^dldd + \ \+ sgn(cr)S" 


VkSte *2 VkSte “2 

H-;;- I lOaOikd' 


-Ed 


— Ek — Ed 


—7o<t7o<t' 
, -£d 


S^^^dldd + (5,^/^ f i + sgn(CT)5'" 


(15) 


Correspondingly, the four parts of the second term of AHll are given by 


(I) 


ill) 


J2 


kk' ,(T(T 

X 


Wfe 


ul' 


Ed 




\vk' 


+ U-Ek^ Ed + U-Ek'J + [Ek+Ed + U^ Ek>+Ed + u) 

1 


S^,^d„dl + S„„, ( - - sgn(a)5' 


- X] 1414 , 

kk' ,(Tij' 


UkVk 


Uk'Vy 


Ed E U — Ek Ek' + £d + t/ 


at at,-, + 

k/T k'cr' 


'^k'^k 


'^k' '^k' 


Ek Sd U Ek' Sd U 


^kcT^k' I 


S„,i^d„dl + S„„, ( - - sgn(a)5'" 


(777) : ^ 14 


UkSte *2 ^ ukSte *2 


+ 


££i + 77 — Ek —Ed — U 
UkSte’’^ 


+ 


I UkSte" ^ 

Y —Ed — U Ed E U — Ek 
'1 


yoo-afco-' 


Ofco- 70 (T' 


VkSte" 


VkSte “2 ^ VkSte *2 


Ek + £d + 77 —Ed — U 


vlSte"-' 


^ku'yOtr' 


-Ed-u"^ Ek + Ed + u] 


Sy^^dadl + S„y ( - - sgn(CT)5" 


i^V): E 




-Ed-U 


7o<t70(t' 


Sy,^dadl + s„y { - - Sgn(cr)5" 


(16) 


Here the following relations have been used, namely, dt^{l — ng)dcr'{l — ngi) = Scr'gdldg + + sgn{a)S^] 

and dynada'ns' = Sa'sd^da + So-'d^ - sgn((T)S'^]. 

In Eqs. (na-diii), we can find that the first parts of them exactly yield the normal Kondo correlation 
between the QD and the s-wave superconductor. By defining the spin density operator S = ^dtjCTaa'da', 
these two parts can be systemized as 


kk' ,<7<7' 


Ai7LZ.i = - E + + 


„ 'y ^ Wkk'oft,ja:k'a + ^ ^ Lkk ') (17) 


where J^-y = VkVyiy^ + 




kk' a 


kk' 




\vk'V 


Ek — Ey —Ed Ek —Ed — U E ^,—£d — U EkESd E}^,-\-£d Ek-^-CdEU E}^,-\-£d~EU 


dl — V, V ,I _ A _ "AL _[iTd_ l"fc'l I l"fcl I \'’k'\ \ Ti 

'Jkk' - ''kVk'(Ek-ed ^ Ey-Sd Ek-Sd-U Ey-Sd-U Ek+ej. Ey+ed ^ Ek+Sd+U ^ Ey+ed+U 

clearly shown that AHll i reflects the normal spin correlation between the localized state in the QD and 
the continuum states in the s-wave superconductor. In the case of negative spin-correlation parameter, the 
two quantum states will correlate with each other in the antiferromagnetic manner, as a result, the Kondo 
effect will come into being. 

The second parts in Eqs. (nnn-disi) describe the singlet-triplet correlations. Their contributions to AHll 
are 

AHll ,2 = - ■*" Afefe')] “*' 9 X/ + h.c., (18) 


kk' 


kk' 


UkVk 


in wbirb A* — V, V,, (—AAk __ uyvy , uyvy n , _ T/ T/, ( i 

m WPlCh - VkVk'KEk-ed-U Ek+Sd Ey+ed + U + Ey-ed> ^ ° ~ ^kVk'{Ek-Sd-U + Ek+ed 

F „ -I- )_ The first term in this equation shows the correlation between the localized state 

in the QD and the Cooper pair in the s-wave superconductor. In the second term, it suggests that the 
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correlation between the localized state in the QD and the Cooper pairs in the s-wave superconductor has an 
opportunity to break up the Cooper pairs and reconstruct the correlation between the localized state and 
new spin-conserved Cooper pairs. Both i and Aid/,/, 2 originate from the coupling between the QD 

and s-wave superconductor, and they have been observed by Salomaaj^ 

Next, the third parts in Eqs. (Ha-dHl) reflect the contribution of the Kondo QD to the Josephson effect. 
It can be found that AHll ,3 = -|- idcross + idr, where 


Hi 


Hr. 


Hj' 


- E + h.c., 

k,(T(T' 

- Y, {MikS.a'S^ + M 2 kCTl,,Sy + + Mik<jy^,S^)ai^joa' + h.c., 

k,(T(T' 

E Tk ,(7(7' +h-c.. 

k^aa' 


(19) 


Relevant parameters here are respectively given by 


Kt = -{Kir 

Mife = iM;^ 
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UkStVki 
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\vk\StVk{ 


Ek-ed Ek- £d-U £d £d + U 

1 111 


-ie /2 


-iei2 


Ek + £d + u Ek + £d £d £d + u 
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with ry-f-j, = and ry^i- = —It clearly shows that due to the presence of the Kondo QD, the 
Josephson effect becomes relatively complicated. The first term, i.e., Hi, exactly describes the Josephson- 
Kondo effect. It can be explained as the Kondo correlation between the localized state and the continuum 
Andreev reflection states, idcross denotes additional spin-exchange correlation. After a further analysis, one 
can find that this term contributes little to the Josephson effect because it only contains the correlations 
between and S'“ ((t“ and S^) {a = x,y). As for Ht, it exactly reflects the screening of the Kondo effect 
induced by Joesphson effect. 

The last terms in Eas. dT5]) - (froi) represent the correlation between the QD and Majorana doublet. Via 
a straightforward derivation, we can simplify them to be Aid /,/,_4 = —JM^na' ■ S'yocrTocr' in which 
Jm = f ^_rr )- It seems that Aid/,/, 4 describes the Kondo-typed correlation between the QD 

and Majorana doublet. However, since Majorana doublet are bound states, such a correlation can only be 
considered to be the normal spin exchange correlation. One can also find that this kind of correlation is very 
weak and it disappears at the position of electron-hole symmetry. In addition, note that due to the special 
anti-commutation relation between the Majorana operators, such a correlation is anisotropic. Up to now 
we have clarified the complicated correlation properties in the case of one Kondo dot embedded in a hybrid 
junction with the Dlll-class TS and s-wave superconductor. 

Based on the above deduction process, it can be clearly found that when the DZ/7-class TS degrades into 
D-class, both AHll ,3 and AHll ,4 will disappear completely. This can be viewed as the essential difference 
between the two kinds of TSs. 


IV. DISCUSSION ABOUT THE JOSEPHSON-KONDO EEEECT 

According to the theory developed in the above section, in the considered hybrid junction, special 
Josephson-Kondo effect can be driven when the QD is in the Kondo regime. In this section, we would 
like to present a comprehensive analysis about such an effect. For convenience, we rewrite Hi as 

Hj^ — ^ ^ ^ J- h.c., 

k,<7<7' 


( 20 ) 
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where /Ci = /C 2 = and /C 3 = k\. Next, we would like to describe the physics picture of the Josephson- 
Kondo effect by defining the Green function 


Gs,w {kr, k'r') = ^ {kr, k'r') = - ^ 


(-ir 



dri • ■ ■ dr„{TT[77i'(Ti) • • • HK{Tn)akr,{T)a\,^,{T')])o. 


The reason is that (/cr, fc'r') reflects respective-order scattering processes which can be described by 

the Feynman diagrams after the Wick contraction. It is easy to find that the first-order scattering process 
is zero because the Wick contraction contains (alg'yos')o and (josaks'}o- Therefore, (/cr,/c'r') makes 

the leading contribution to the Josephson-Kondo physics with 






(n)7os'i (ri)4i (n) 


^ 11^2 sis^,S 2 S^ 


+^i(o-s;si ■o-!.i,.i)> 0 s;(n)afcisi(n)d^,(ri)d,,i(Ti)] x 42S2 4^4os'J'^2)dtjT2)d^^(T2) 

+^i(o- 4 s 2 ■^^i^2yyos'J'^2)ak2S2(r2)dl^(T2)d^^(T2)] X afc„(r)at,^,(r')})o. 

With the help of Wick theorem, we plot the connected Feynman diagrams which reflect the second-order 
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FIG. 2: Feynman diagrams of the second-order scattering process. 


scattering process in Fig.2. Here the solid, dashed, and dotted lines correspond to the free propagators 
G^slrjiik,uJn), and (w/), respectively. In comparison with these two diagrams, one can 

be sure that the scattering process involved in the first one contribute dominantly to the Josephson-Kondo 
effect. 

In the following, we concentrate on the calculation about the contribution of the first Feynman diagram 
in Fig.2 to the Josephson-Kondo effect. It should be noted that the calculation about the contribution of 
the Feynman diagram is to evaluate its-induced self-energy. According to the Dyson series in random-phase- 
approximation theory, the total particle-particle interactional vertex can be given by 




/32 

kO) 


p ^ iOs^UJi S2 ,^' 2 ; S2,^'2 


xF° 

S2l'2S'2^'2'~' M,82'^ 




where = \lCi{cr accompanied by its illustration in Fig.3. 

Next, suppose 

= F 3sis[4iI2[ + 'y ^ F i{(^sis[ ■ O’!/!!/') > 


( 21 ) 


( 22 ) 


we can reexpress the above Dyson equation, i.e., 


+ cos9+ hic4^r]+ \\F\^ ( 23 ) 



I 



FIG. 3: Illustration of the Dyson series of the vertex in the random-phase approximation. 


in which j G^^\uJi)G^^\uJm-i^i)- Related deduction 

process can be referred in Appendix A. Substitute Ea. d^ into Ea. d^ . we get the following result 


r° = (^|/c-^pcos 0 r° + f 

i 

+ + I|/C^pcos6irf)n‘5^(w„)n^^(a;„) 

with r'a = Fp + F}y (e ap-y = 1). Although four variables exist in these two equations, by summing up F] 
we can solve these equations and figure out F°, i.e., 


F°{0Jm) = 


Y:iFi 


[1 - (||K:-L|2 + i|/C-L|2cos6i)nsen*fe][i _ ll/c^pcosenS-PnA^e] _ ||/C-L|4cos6»(nsenA^e)2- 


According to this result, the expression of F} can be clarified. It is evident that the solution of F°, F] and 
then r depends on the calculation of (See Appendix B). 

With the help of the above discussion results, the conductance quasiparticle irreducible self-energy can be 
calculated by using fluctuation-dissipation theorem. As a consequence, it obeys the following equation 


SK) = ^E' 




/ + 00 1 ^ ^ 1 

duj{ -Im[n^^(u;)r(a;)]}- E T- -:-^ 

-oo TT P ^ {lUJn - lUJmj 


){iuJm - Uj) ’ 


(24) 


in which r 6,^,^ = 




■jyO-s',s 2 After 


the Matsubara-frequency summation on ujm (w™ = 27rm//3 and m G the relaxation time that resembles 
the Josephson tunneling conductance can be estimated, since it is related to the imaginary part of self-energy 
E with Tp{ 9 ) ^ —ImE(a; + 


V. SUMMARY 

In summary, we have presented the discussion about the spin-correlation effect in the QD-embedded 
Josephson junction formed by the indirect coupling between one Dlll-class TS and s-wave superconductors. 
By carrying out the Schrieffer-Wolff transformation, we have found that the U/JJ-class TS contributes 
additional anisotropic spin correlations to the modification of Kondo physics, which are related to the 
Josephson phase difference between the superconductors. These special spin correlations originate from 
the continuous Andreev reflection between the superconductors provides a new continuum states for such 
a kind of spin correlation. What’s notable is that if the time-reversal invariance in the TS is broken, such 
terms will disappear. These results reflect the fundamental difference between the Majorana bound states 
in the D- and Dlll-dass TSs. As a typical case, the Josephson-Kondo effect has been analyzed with the 
perturbation method, and relevant parameters have been calculated. We believe that this work is helpful 
for understanding the Kondo effect in the system with Dill-class TS. 
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Appendix A 

In the solution procedure of F, the following step should be involved, i.e., 


1 1^1 
- = tE - t(/c^/cii* + • a,.,,. 




which iEa"/ 3 "( 


’■/3/3")( 




^r- I -=2 = i-'^2‘'i 2'"' ' (25) 

/3"/3') = jSaa'Si3i3> — ^cTaa' ' O'/?/?'- Higher-Order spin summation 


can be derived in the same way. It is 

g ^ ^ ) [ICj {o’sys!, ' ^1221^2') ][]Ck{<y siS 2 ' )1 


S2l^2 s'ryU'ry 


lS2 1-^2/ 


= E[||^^I^^'>2s'/-'2U - cosie + ■ (T^^^,^][K:k{cTs^^ 

^21^2 

= ~ ;jl^^l^cos( 6 i + 7 r) E K.k{.cTs^s2 ■ 0-1^1122)'° {a-S2S[ ■ 

S2i^2 

= cos{e + n)'^JCkSs^s{S^i,.{ - cos{e + n)K:'k{(T^^„>^ •o'.-iU)* 


3 1 

~ '7^0^3132^1^112'^ ~ '^ 12112 '^) , 


(26) 


^-.J-0102 2 

where Kq = /Ci) cos 0 and = — ||/C^p/Ci — l|/C^p/C'cos 0 with IC'^ = IC 2 = IC^ + and 

/C^ = IC^). 


Appendix B 


The pair-bubble diagram corresponds to the Matsubara-frequency summation, and it can be calculated 
as follows 

rr-SP/ \r<P( ^ fi^p) 

H (oJm) — a / . ^0 (Pi ^sJGq (uJm l^.s) — . + . p ^ ’ 

P Ep liOm Ep Sd 


n“^(.;^) = iE^o""(^i)Go"’(a;™-^*) = 4 E— Y-^- ■ 

P ^ P ^ lijJi t{uJm - Wj) - Sd 


(27) 


By performing the analytical continuation, they can be transformed into 


n-(.) = EE 


1 - i tanh (ided) 1tanh (Id-Fp) e^=o p{E)dE 1 1+^ p(Fi) tanh (1/3-F) 




iu) — (£/p + £d) 2 Ep £ci ~ 


ioj — E 2 


E — iui 


-dE 


»;+i0+ 




r+D 


-D 


Pi^BE . 2 ,1^ f+^ p{E)tanh{^pE) in 


^-E -*-p(^) + E 


'-D 


„ dE -\ -tanh 

E — uj 2 


^)p(a^), 


r^MPl N f{iu!m—£d) /(O) Sd=0 1 ii^-»i^+iO+ 1 

ii [LO) — -:-h ^-/■ ^-/■ - .„ I . 

^ H” ^0 


(28) 


due to p{E) = I'^l) and v = sgn(if -|- lAI). In this derivation process, we set Ed to zero in order to 

V-® 1^1 

eliminate those non-physical states in pseudofermion representation. 
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